AD  A  131179 


zty.e  jot /88 


DNA  5991 F 

EARLY  TIME  STARFISH  CALCULATIONS 


F.  E.  Fojen 
R.  W.  Kilb 


Mission  Research  Corporation 

735  State  Street 

Santa  Barbara,  California  93101 


1  February  1982 


Final  Report  for  Period  30  April  1979—15  June  1980 


CONTRACT  No.  DNA  001-79-C-0276 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


DTIC 

SELECTE 

AUG  8  1983 

<£  B 


THIS  WORK  WAS  SPONSORED  BY  THE  DEFENSE  NUCLEAR  AGENCY 
UNDER  RDT&E  RMSS  CODE  B322079464  S99QAXHC06243  H2590D. 


Prepared  for 
Director 

DEFENSE  NUCLEAR  AGENCY 
Washington,  DC  20305 


88  OG 


(04 


UNCLASSIFIED 


SECuRtTy  Cl  ASSJPlC  ATJON  OF  T  mi  PAGE  '»hm  Dale  Fmered, 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

HEEORE  r~OMPEHT!NCi  FORM 

1  «EPjRTNjW'tU  :  GOVT  ACCESSION  NO. 

DNA  5991F  /)  /  /  A 

|  3  RECiPlENT-t  C  AT  Al.'I'j  NUMBER 

4  T|TlF  and  Subtitle'. 

i  Type  Of  REPORT  &  PERIOD  covered 

Final  Report  for  Period 

30  Apr  79  to  15  Jun  80 

EARLY  TIME  STARFISH  CALCULATIONS 

6  PERFORMING  OR3  REPORT  NUMBER 

MRC-R-572 

A  .  'r  -. 

F.  E.  Fajen 

R.  W.  Ki lb 

8  contract  or  grant  NUMBER*' 

DNA  001 -79-C-0276 

•)  Pf  R  c  R  M '  N  C  C'R  j  AN;  J  A  T  -  CN  *(«M;  A  N  G  ADDRESS 

Mission  Research  Corporation 

735  State  Street 

Santa  Barbara,  CA  93101 

PROGRAM  E_EMEn~  PROJECT  A s*f 

AREA  6  JMT  NUMBERS 

Subtask  S99QAXHC062-43 

1*  *C\*R.\UPI3  O  F  F  l  "  f  NAME  and  ADDRESS 

1 2  BEP3»T  I«TE 

Di rector 

1  February  1982 

Defense  Nuclear  Agency 

Washington,  D.C.  20305 

1  3  NUMBER  O  F  PAGES 

80 

'f  0.  S-Ri9  .  T ‘  ON  STATEMENT 


Approve*,  for  public  release;  distribution  unlimited. 


***  St4’FMF.Nt  'of  »/:«■  flhs 


10  SuPPlEMEN’ARv  notes 


This  work  was  sponsored  by  the  Defense  Nuclear  Agency  under  RDT&E  RMSS 
Code  B32 2079464  S99QAXHC06243  H2590D. 


19  KEv  WORDS  'Continue  <>r?  reverse  side  if  necessary  and  identify  bv  block  number i 

Nuclear  Burst  above  150  km.  Debris  Energy  Patches 

Expansion  from  0  to  1  sec.  Conjugate  Patches 

Blast  Wave  Behavior  Velocity  Spectrum  of  Conjugate 

Collisionless  Magnetohydrodynamics  Patches 


20  ABSTRACT  'Continue  on  reverse  side  if  necessary  and  identify  h\ •  block  number) 

In  this  report,  we  present  results  from  recent  CMHD  (Collisionless 
Magnetohydrodynamics)  code  calculations  of  the  early  time  [  0  to  1  sec.) 
evolution  of  the  Starfish  event.  The  emphasis  of  this  report  is  on  the 
velocity  spectrum  of  the  energy  going  to  the  conjugate  patches.  We  present 
a  parameter  study,  utilizing  results  from  eight  simplified,  "CMHD-like" 
calculations,  which  demonstrates  the  dominant  role  played  by  ionizing 
collisions  during  the  blast  wave  expansion. 


DD  1473  EDITION  OF  1  NOV  *5  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  f^e n  Hat*  Entered 


Sect  ion 


TABLE  OF  CONTENTS 


LIST  OF  ILLUSTRATIONS 

1  INTRODUCTION 

2  THE  CMHD  CODE 

3  RECENT  CMHD  CALCULATIONS  OF  STARFISH 

4  PATCH  SPECTRUM  SENSITIVITY  STUDY 

REFERENCES 


LIST  OF  ILLUSTRATIONS 


Figure 


The  CMHD  computational  grid  used  for  Starfish.  Radial  11 
lines  emanating  from  the  burst  point  are  Eulerian 
surfaces,  while  the  vertical  lines  which  represent  mag¬ 
netic  flux  surfaces,  are  Lagrangian. 

The  CMHD  computational  grid  at  .025  sec.,  .03  sec.,  13 

.05  sec.,  and  .10  sec.,  during  a  simulation  of  the 
Starfish  event. 

The  CMHD  computat ional  grid  at  .20  sec.,  .40  sec.,  .60  14 

sec.,  and  .80  sec.,  during  a  simulation  of  the  Starfish 
event . 

Comparison  of  Starfish  NCR  power,  as  inferred  from  19 

optical  data,  with  f(t)  =  t'1 . 

Comparison  of  Ne  contours  from  the  three  Starfish  24 

simulations.  The  values  contoured  are  in  units  of 
electrons/cm3.  The  contour  lines  terminate  at  the 
surface  of  the  burst-produced  magnetic  bubble  and  at 
the  edge  of  the  computational  grid. 

Comparison  of  N„  contours  from  the  three  Starfish  25 

simulations.  The  values  contoured  are  in  units  of 
electrons/cm  .  The  contour  lines  terminate  at  the 
surface  of  the  burst-produced  magnetic  bubble  and  at 
the  edge  of  the  computational  qrid. 

Energy  delivered  to  the  downward  (NCR)  debris  energy  26 
patch  as  a  function  of  cross-field  radius. 

Comparisons  of  fast  ion  contribution  to  downward  (NCR)  29 
debris  energy  patch  spectra. 

Comparison  of  fast  neutral  contribution  to  downward  30 

(NCR)  debris  energy  patch  spectra. 

Comparison  of  downward  (NCR)  debris  energy  patch  31 

spectra. 


2 


LIST  OF  ILLUSTRATIONS  (Continued) 


Figure  Page 

11  Geometry  of  flux  tube  no.  5  at  various  times  during  33 

the  blast  wave  evolution. 

12  Energy  and  mass  sources  for  the  entire  flux  tube,  as  34 

determined  from  the  full  CMHD  run. 

13  Time  integrated  values  of  the  energy  and  mass  sources.  35 


14  Schematic  representat ion  of  the  flow  of  energy  and  mass  37 
in  the  single  flux  table  calculations.  Angular  sectors 

5  through  21  are  "active"  for  this  flux  tube  -  sectors 
1  and  22  are  straight  down  and  straight  up,  respect i vely, 
where  the  non-orthogonality  of  the  CMHD  coordinate  sys¬ 
tem  causes  singularities,  and  sectors  2  through  4  inter¬ 
sect  this  particular  flux  tube  below  the  bottom  of  the 
CMHD  grid.  S  represents  the  action  of  the  SCATTER 
subroutine,  which  takes  energy  from  transverse  thermal 
and  puts  it  into  parallel  streaming,  and  T  represents 
the  transport  of  parallel  streaming  energy  up  and  down 
the  flux  tube. 

15  An  example  of  the  initial  ion  distribution  function,  39 

dM 

-gy^j-,  within  one  angular  sector. 

dM 

16  An  example  of  the  distribution  function,  ^gy,  within  40 

one  angular  sector,  after  the  diffusive  SCATTER  sub¬ 
routine  has  deflected  energy  from  the  transverse  direc¬ 
tion  into  the  parallel  direction. 

dM 

17  An  example  of  the  distribution  function,  -ggy,  after  the  41 

WP  SCATTER  routine  has  deflected  energy  from  the  trans¬ 
verse  direction  to  the  para1 lei  direction. 

18  Typical  contours  of  the  logarithm  of  the  distribution  42 

of  ion  mass  per  unit  parallel  velocity,  with  dM/dUn  in 
units  of  gm/(cm/sec),  as  a  function  of  parallel  velocity 
and  angular  sector  number.  The  outermost  contour  has 
dM/dUH  =  10*7  gm/(cm/sec),  and  the  innermost  contour 

has  dM/dUn  =  10'3  gm/(cm/sec)  in  this  example  at  a  time 
0.1  sec. 


3 


LIST  OF  ILLUSTRATIONS  (Continued) 


Figure  Page 

19  Transverse  energy,  parallel  energy,  and  patch  energy  44 

as  functions  of  time  for  CASE  1. 

20  Distribution  function  contours  for  CASE  1  at  times  of  46 

.026,  .030,  .050,  and  .100  seconds. 

21  Distribution  function  contours  for  CASE  1  at  .20,  47 

.40,  .60,  and  1.00  seconds. 

dM 

22  The  velocity  spectrum,  -wr-,  for  the  downward  patch  49 

for  CASE  1.  au" 

23  Transverse  energy,  parallel  energy,  and  patch  energy  50 

as  functions  of  time  for  CASE  2. 

24  Distribution  function  contours  for  CASE  2  at  times  of  51 

.026,  .030,  .050,  and  .100  seconds. 

25  Distribution  function  contours  for  CASE  2  at  times  of  52 

.20,  .40,  .60,  and  1.00  seconds. 


26  The  velocity  spectrum,  -ttt-,  for  the  downward  patch  53 

for  CASE  2.  dU" 

27  Transverse  energy  as  a  function  of  time  for  CASE  2,  55 

CASE  3,  CASE  4,  and  CASE  5. 

28  Parallel  energy  as  a  function  of  time  for  CASE  2,  57 

CASE  3,  CASE  4,  and  CASE  5. 

29  Patch  energy  as  a  function  of  time  for  CASE  2,  CASE  3,  58 

CASE  4,  and  CASE  5. 


30  Comparison  of  distribution  function  contours  for  59 

CASE  2,  CASE  3,  CASE  4,  and  CASE  5  at  .05  seconds. 

31  Comparison  of  distribution  function  contours  for  60 

CASE  2,  CASE  3,  CASE  4,  and  CASE  5  at  .40  seconds. 

32  Comparison  of  distribution  function  contours  for  61 

CASE  2,  CASE  3,  CASE  4,  and  CASE  5  at  1 .0  seconds. 

33  Comparison  of  patch  spectra  for  CASE  2,  CASE  3,  63 

CASE  4,  and  CASE  5. 


4 


LIST  OF  ILLUSTRATIONS  (Concluded) 


Figure 

Page 

34 

Transverse  energy  as  a  function  of  time  for  CASE  2, 
CASE  6,  CASE  7,  and  CASE  8. 

64 

35 

Parallel  energy  as  a  function  of  time  for  CASE  2, 

CASE  6,  CASE  7,  and  CASE  8. 

65 

36 

Patch  energy  as  a  function  of  time  for  CASE  2, 

CASE  6,  CASE  7,  and  CASE  8. 

66 

37 

Comparison  of  distribution  function  contours  for 

CASE  2,  CASE  6,  CASE  7,  and  CASE  8  at  .05  seconds. 

67 

38 

Comparison  of  distribution  function  contours  for 

CASE  2,  CASE  6,  CASE  7,  and  CASE  8  at  .40  seconds. 

68 

39 

Comparison  of  distribution  function  contours  for 

CASE  2,  CASE  6,  CASE  7,  and  CASE  8  at  1.00  seconds. 

69 

40 

Comparison  of  patch  spectra  for  CASE  2,  CASE  6, 

CASE  7,  and  CASE  8. 

71 

5 


SECTION  I 
INTRODUCTION 

For  nuclear  bursts  occurring  above  about  150  km  altitude  and 
below  about  600  km,  an  important  aspect  of  the  burst  phenomenology  is  the 
debr'is  energy  patch.  Analysis  of  photographic  data  taken  following  the 
Starfish  nuclear  test  has  indicated  that  one-fourth  to  one-half  of  the 
kinetic  yield  of  the  weapon  arrived  in  each  magnetic  conjugate  region 
during  the  first  few  seconds  after  the  burst. 

The  existence  of  the  debris  energy  patch  was  predicted  prior  to 
the  Starfish  event  by  longmire,  et.  al.,1  and  their  prediction  proved  to 
be  essentially  correct.  The  current  understanding  of  the  phenomenon  is 
that  weapon  debris  ions  and  air  ions  that  are  accelerated  to  high  veloc¬ 
ities  (107  cm/sec  to  a  few  times  108  cm/sec)  by  the  expanding  debris-air 
blast  wave,  stream  parallel  to  the  (distorted)  geomagnetic  field  lines 
from  the  burst  region  to  each  magnetic  conjugate,  where  they  are  stopped 
by  ordinary  atomic  collisions  and  thereby  deliver  their  energy  to  the 
background  air,  heating  and  ionizing  it.  Ionization  thus  produced  in  the 
debris  energy  patches  may  persist  for  several  hours  after  the  burst.  It 
is  likely  to  develop  structure  -  i.e.,  striations  -  and  it  therefore  is 
expected  to  degrade  the  performance  of  communication  links  which  propagate 
through  it. 

Under  DNA  funding,  a  highly  specialized  computer  code  was  devel¬ 
oped  at  Mission  Research  Corporation  (MRC)  to  treat  the  early  time  phenom¬ 
enology  of  nuclear  bursts  in  the  150  km  to  600  km  altitude  range.  This 
code,  known  as  CMHD  (Col 1 isionless  Magnetohydrodynamics) ,  briefly 
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described  in  Section  2  of  this  report,  is  fully  described  and  documented 
in  References  2  and  3. 


Seven  high  altitude  burst  calculations,  covering  a  range  of 
burst  altitudes  and  yields,  have  been  completed  with  CMHD,  and  have  been 
published  in  reference  3.  The  first  of  these  was  a  simulation  of 
Starfish.  Comparisons  of  predicted  northern  conjugate  region  (NCR)  power 
as  a  function  of  time  after  burst  were  made  with  the  published  power-time 
curve  inferred  from  the  photographic  data,  and  adequate  agreement  was 
obtained. 


However,  the  Starfish  photographic  data  do  not  yield  sufficient 
information  to  verify  all  aspects  of  the  CMHD  code.  And,  unfortunately, 
the  KLYSMA  -  SCORPIO  -  CYLWHAT  family  of  codes  at  the  Naval  Research  Lab¬ 
oratory  (NRL),  when  applied  to  Starfish,  have  given  results  which  differ 
in  certain  respects  from  those  of  CMHD.  In  particular,  the  two  sets  of 
results  differ  with  regard  to  the  velocity  spectrum  of  the  fast  ions  which 
stream  from  the  burst  area  to  the  conjugate  regions  to  generate  the  debris 
energy  patches.  The  NRL  spectrum  tends  to  be  harder,  to  cause  the  patch 
ionization  to  occur  at  lower  altitudes,  where  molecular  deionization  can 
proceed  rapidly,  and  thus  produce  less  persistent  ionization  than  does  the 
CMHD  spectrin  (if  the  energy/cm2  of  the  deposition  along  the  field  lines 
is  less  than  -107  ergs/cm2  so  that  molecular  species  still  predominate 
after  the  NRL-type  energy  deposition). 

This  report  presents  results  of  two  additional  Starfish  runs 
that  have  been  made  with  the  CMHD  code  as  part  of  the  continuing  effort  to 
resolve  this  disagreement  between  MR C  and  NRL  with  regard  to  the  debris 
energy  patch  spectrum.  These  additional  CMHD  calculations  are  described 
in  Section  3. 
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In  Section  4,  several  simplified  "CMHO-like"  calculations  are 
presented,  which  demonstrate  the  degree  to  which  the  debris  energy  patch 
spectrum  is  sensitive  to  different  aspects  of  the  physics  in  the  CMHD 
code.  The  result  of  this  study  is  that  the  spectrum  is  highly  sensitive 
to  the  rate  at  which  background  air  is  ionized  by  collisional  processes  in 
the  expanding  blast  wave.  The  detailed  treatment  given  to  these  colli¬ 
sional  processes  in  the  CMHD  code,  and  the  lack  of  a  corresponding 
treatment  in  the  NRL  codes  is  believed  to  be  the  primary  source  of  the 
discrepancy  between  the  MRC  and  NRL  results. 
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SECTION  2 
THE  CMHD  CODE 


The  physics  embodied  in  the  CMHD  code,  the  derivation  of  the 
mathematical  equations,  and  the  code  itself  are  documented  in  references  2 
and  3.  The  CMHD  code  is  believed  (at  MRC)  to  represent  the  current  state 
of  the  art  for  early  time  nuclear  burst  phenomenology  in  the  relevant 
altitude  regime.  A  very  brief  synopsis  of  the  CMHD  physics  now  follows. 

Ideally,  one  would  like  to  solve  the  Boltzmann  transport  equa¬ 
tion  for  each  ion  and  neutral  species,  resulting  in  several  distribution 
functions,  each  a  function  of  time  and  six  space  and  velocity  variables. 
However,  this  is  well  beyond  the  capacity  of  present  day  computers,  and  we 
are  forced  to  eliminate  some  of  the  independent  variables.  One  space 
variable  and  two  velocity  variables  can  be  eliminated  by  assuming  azimuth¬ 
al  symmetry  -  appropriate  for  the  case  of  an  initially  vertical  geomag¬ 
netic  field  -  and  by  assuming  that  all  ions  in  the  same  spatial  location 
have  the  same  average  velocity  in  the  direction  perpendicular  to  B.  The 
latter  assumption  is  justified  on  the  basis  of  momentum  coupling  in  the 
transverse-to-^  direction,  which  is  generally  believed  to  occur  on  a 
length  scale  of  an  ion  gyro  radius  or  shorter.  Appropriate  moments  of  the 
Boltzmann  equation  are  used,  then,  to  derive  a  set  of  equations  which  are 
fluid-like  in  the  transverse-to-^  direction,  but  which  retain  the  arbi¬ 
trary  distribution  function  in  parallel  velocity  space.  That  is,  ions 
move  together  as  a  fluid  in  the  perpendicular  direction,  but  are  allowed 
to  stream  independently  in  the  parallel  direction. 
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Altitude  (km) 


The  CMHD  computational  grid  used  for  Starfish.  Radial  lines 
emanating  from  the  burst  point  are  Euler ian  surfaces,  while 
the  vertical  lines  which  represent  magnetic  flux  surfaces,  are 
Lagrangian. 
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Figure  1. 


The  CMHD  computational  grid,  illustrated  in  Figure  1,  consists 
of  cells,  enclosed  by  two  sets  of  surfaces.  The  radial  surfaces  emanating 
from  the  burst  point  remain  fixed  in  space  (Eulerian  surfaces)  as  the  ions 
move  upward  or  downward  through  them  from  one  cell  into  a  neighboring 
one.  The  other  set  of  surfaces,  initially  vertical,  represent  magnetic 
flux  surfaces,  and  move  outward  with  the  ions  (Lagrangian  surfaces)  as 
they  expand  away  from  the  immediate  vicinity  of  the  burst.  Initially, 
each  cell  contains  ions  present  in  the  ambient  ionosphere  plus  ioi.s 
created  by  X-rays  from  the  burst.  Further  ionization  may  be  generated  by 
deposition  of  UV  energy  and  by  collisional  ionization  by  ions  or  electrons 
of  background  neutrals.  This  collisional  ionization  of  neutrals  is  an 
important  source  of  ions,  which  can  cause  an  exponential  increase  in  the 
ion  mass  within  a  flux  tube  as  the  blast  wave  expansion  progresses. 

Figures  2  and  3  show  the  configuration  of  the  CMHD  grid  at 
several  times  during  the  blast  wave  expansion  for  one  of  our  recent  Star¬ 
fish  simulations.  At  .025  seconds,  the  problem  initialization  time,  only 
the  two  innermost  flux  surfaces  are  distorted.  As  the  expansion  proceeds, 
the  volume  of  the  magnetic  bubble  increases  dramatically,  and  by  0.4 
seconds,  all  of  the  flux  surfaces  are  distorted,  some  quite  severely. 
Simulations  such  as  the  one  depicted  here  are  usually  run  out  to  a  problem 
time  of  one  to  two  seconds  after  burst  time. 

As  a  cell  is  compressed  and  set  in  motion,  the  ions  within  that 
cell  acquire  "transverse  thermal"  energy  as  they  are  picked  up  by  the  mov¬ 
ing  field  lines  and  gyrate  about  these  field  lines.  There  may  also  be 
plasma  turbulent  processes  which  excite  transverse  thermal  modes,  but  such 
processes  are  not  treated  in  detail  by  CMHD.  Rather,  the  code  solves  the 
fluid-like  equation  derived  from  the  appropriate  moment  of  the  Boltzmann 
equation  to  obtain  the  ion  transverse  thermal  energy.  If  the  resulting 
transverse  pressure  is  sufficiently  large  compared  to  the  parallel  pres¬ 
sure,  then  the  plasma  is  unstable,  and  one  or  more  plasma  instabilities 
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Figure  2.  The  CMHD  computational  grid  at  .025  sec.,  .03  sec.,  .05  sec 
and  .10  sec.,  during  a  simulation  of  the  Starfish  event. 


The  CMHD  computational  grid  at  .20  sec.,  .40  sec.,  .60 
and  .80  sec.,  during  a  simulation  of  the  Starfish  event 


will  turn  on,  deflecting  ion  thermal  velocities  from  the  transverse  direc¬ 
tion  into  the  parallel  direction.  The  turn-on  condition  for  the  mirror 
instability4  is 
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In  the  previously  published  CMHD  calculation,  this  instability  was  simula¬ 
ted  by  the  SCATTER  subroutine,  by  performing  velocity  space  diffusion  to 
increase  P  at  the  expense  of  P  whenever  the  above  condition  was  met. 
The  result  of  this  was  to  drive  the  velocity  distribution  toward  a  Maxwel¬ 
lian,  which  is  generally  what  we  would  expect  from  turbulent  processes. 
Other  forms  of  the  SCATTER  subroutine,  simulating  other  types  of  scatter¬ 
ing  in  velocity  space,  have  been  exercised  in  recent  CMHD  calculations,  as 
will  be  described  in  Sections  3  and  4. 


As  the  ions  follow  their  independent  paral lei -to-8  trajectories, 
they  undergo  several  interactions.  Ions  with  different  parallel  veloci¬ 
ties  interact  via  ordinary  coulomb  collisions.  This  is  generally  impor¬ 
tant  only  at  relatively  large  ion  densities  such  as  when  the  bottom  of  the 
grid  is  below  200  km  altitude,  and  ions  created  at  higher  altitudes  find 
themselves  streaming  through  essentially  stationary  ions  created  by  X-rays 
and  UV.  Ions  interact  with  the  magnetic  field,  via  the  mirror  force, 
which  accelerates  them  from  regions  of  high  field  strength  toward  regions 
of  low  B,  consituting  an  interchange  between  transverse  thermal  energy  and 
parallel  streaming  energy.  As  ions  stream  through  a  region  of  curved 
field  lines  which  is  in  motion,  there  is  interchange  between  transverse 
streaming  energy  and  parallel  streaming  energy  (sling  effect).  Charge 
exchange  between  a  moving  ion  and  a  stationary  neutral  can  occur,  result¬ 
ing  in  a  stationary  ion  (which  is  then  swept  outward  along  with  the  moving 
magnetic  flux  surface)  and  a  fast  neutral.  The  resulting  fast  (or  "hot") 
neutrals  are  then  treated  with  transport  equations  similar  in  complexity 


to  the  ion  equations.  The  neutrals  may  become  col  1 isional ly  re-ionized, 
or  they  may  stream  out  of  the  computational  grid  as  neutrals.  The  neu¬ 
trals  which  exit  the  grid  in  the  downward  direction  constitute  an  impor¬ 
tant  source  of  deposition  energy  in  the  patch  region  below  the  burst.  The 
fast  neutrals  which  exit  upward  from  the  grid  have  more  than  sufficient 
energy  to  escape  into  interplanetary  space,  thereby  reducing  the  effective 
kinetic  yield  of  the  burst. 

Electrons  are  heated  through  coulomb  collisions  with  ions,  and 
also  by  plasma  turbulent  processes,  to  temperatures  that  would  approach 
several  keV,  if  there  were  not  competing  processes  which  cool  them.  These 
cooling  processes  are  (i)  ionization  of  background  neutrals,  mentioned 
above  as  an  important  source  of  ion  mass,  (ii)  further  ionization  of  air 
ions,  which  can  result  in  significant  populations  of  multiply  ionized 
states  of  oxygen  and  nitrogen,  and  { i  i  i )  electronic  excitation  of  these 
air  ions,  followed  by  radiation  of  ultraviolet  (UV)  photons.  At  altitudes 
below  300  km,  the  UV  radiation  is  a  significant  mechanism  for  cooling  the 
blast  wave,  as  well  as  for  producing  ionization  in  regions  outside  the 
current  position  of  the  blast  wave. 

To  summarize,  then,  the  CMHD  code  was  designed  and  developed  to 
compute,  for  high  altitude  bursts: 

(i)  blast  wave  expansion, 

(ii)  transport  of  energy  away  from  the  blast  wave  by  fast  ions 
streaming  parallel  to  tl,  and  the  various  interactions 
between  these  streaming  ions  and  the  blast  wave, 

(iii)  charge  exchange  losses,  and  subsequent  transport  and  reion¬ 
ization  of  the  resulting  fast  neutrals,  and 

(iv)  multiple  ionization  state  chemistry,  the  associated  product¬ 
ion  of  UV  radiation,  and  the  subsequent  deposition  of  this 
UV  energy. 
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These  processes  are  highly  intertwined  and  interdependent.  Experience  has 
shown  that  it  is  necessary  to  treat  them  in  a  self-consistent  manner  in 
order  to  obtain  a  reliable  computation  of  the  energy  transport  and  of  the 
spatial  location  of  final  energy  deposition  for  a  high  altitude  nuclear 
burst,  which  is  the  objective  of  the  CMHD  code. 
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SECTION  3 

RECENT  CMHD  CALCULATIONS  OF  STARFISH 


The  first  full-scale  nuclear  burst  simulation  completed  with  the 
CMHD  code  was  the  Starfish  calculation  reported  in  reference  3.  This 
calculation,  in  common  with  the  other  six  burst  simulations  that  were 
presented  in  the  same  technical  report,  gave  a  velocity  spectrum  for  the 
debris  energy  patch  which  had  the  approximate*  functional  form 


where  is  the  energy  per  unit  velocity  delivered  to  the  patch  by  ions 
with  parallel  velocity  between  Un  and  Un  +  dll .  This  general  form 
is  consistent  with  the  time  dependence  of  the  observed  Starfish  total  NCR 
power,5  which  is  plotted  in  Figure  4.  Also  plotted  in  Figure  4  is  a  curve 
representing  the  function  f(t)  =  t~l  .  It  is  seen  that  the  two  curves 
agree  quite  well,  for  times  after  the  peak  in  the  power-time  curve. 
Assuming  that  the  time  of  arrival  of  the  fast  ions  in  the  patch  region  is 
inversely  proportional  to  the  velocity,  we  have 


and  thus 


(3-1) 

(3-2) 


*  In  reference  3,  the  patch  spectrum  is  fit  with  a  much  more  elaborate 
function,  which  takes  into  account  the  spatial  dependence  of  the 
spectrum.  The  present  discussion  is  primarily  qualitative,  so  the  more 
approximate  functional  form  used  here  is  quite  adequate. 
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Figure  4.  Comparison  of  Starfish  NCR  power,  as  inferred  from  optical 
data,  with  f( t)  =  t  . 
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which  can  be  combined  with  the  observed  time  dependence  of  the  NCR  power. 


dE 

cTF 


(3-3) 


to  yield  an  estimate  of  the  form  of  the  velocity  spectrum,  as  follows: 


or, 


dE 

an 

dE 

m 


dE  dt  i 

nr  *  au  1 


,-2 


(u-1  )-l»ir2  =  u~l 


(3-4) 

(3-5) 


Thus,  it  is  seen  that  the  Starfish  data  are  consistent  with  a  spatially 
integrated  spectrum  which  is  approx imately  the  same  as  that  predicted  by 
CMHD.  Although  the  exact  computed  spectral  shape  is  a  slowly  varying 
function  of  spatial  location,  the  general  form  of  equation  3-5  is  approxi¬ 
mately  applicable  at  all  radii. 


The  NRL  results,  on  the  other  hand,  show  considerable  dependence 
on  cross-field  radius.  At  small  radii,  the  NRL  computed  spectrum  is  more 
like  -  U+1  .  On  a  spatially  integrated  basis,  however,  their  results 
have  also  been  shown  to  be  approximately  consistent  with  the  optical 
data.  Kilb3  compared  the  NRL  and  MRC  results  with  each  other  and  with 
data  in  considerable  detail,  and  concluded  that  in  spite  of  the  fact  that 
the  velocity  spectra  predictions  by  the  two  groups  are  quite  different, 
and  even  though  the  data  seem  to  imply  that  the  form  of  equation  3-5  is 
approximately  correct,  one  could  not  make  an  iron-clad  distinction  between 
the  different  results  on  the  basis  of  the  published  optical  power  data, 
which  is  spatially  integrated  over  the  NCR  region. 


It  has  therefore  been  necessary  to  attempt  to  resolve  the  issue 
on  theoretical  grounds,  which  has  proved  to  be  a  difficult  task.  As  part 
of  this  on-going  effort  to  resolve  these  discrepancies  between  the  NRL  and 
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MRC  code  results,  we  have  made  two  additional  CMHD  runs  for  the  Starfish 
event.  We  shall  refer  to  the  original  calculation  as  “run  A,"  and  to  the 
more  recent  calculations  as  "run  B"  and  "run  C."  Run  A  used  a  form  of  the 
SCATTER  subroutine  which  performed  velocity  space  diffusion  to  drive  the 
local  velocity  distribution  toward  isotropization  whenever  the  conditions 
for  the  mirror  instability  were  met.  Dr.  Robert  Clark,  Dr.  Peter 
Palmadesso,  and  other  NRL  scientists  objected  strenuously  to  this 
procedure.  The  NRL  philosophy  concerning  the  scattering  of  energy  from 
transverse  thermal  modes  into  parallel  stream  modes  is  expressed  by  the 
following  quotation:6 

"We  focus  our  attention  on  the  coupling  region  for 
the  present.  As  we  have  said  it  contains  a  highly 
compressed  magnetic  field  and  a  mixture  of  hot  air  and 
debris.  The  debris  ions  have  streaming  energies  predom¬ 
inantly  perpendicular  to  the  field.  Coupling  instabil¬ 
ities  in  this  medium  convert  this  energy  into  random 
transverse  motion,  heating  the  debris  and  air  ions 
anisotropical ly  in  the  transverse  direction  so  that 
T^  /T  ~  3.  In  this  way  a  large  part  of  the  transverse 
streaming  energy  of  the  debris  ions  is  transferred 
quickly  to  the  parallel  direction  and  shared  with  the 
ambient,  ions. 

This  process  results  in  a  peculiar  velocity  distri¬ 
bution,  deficient  in  specific  high  velocity  values  at 
which  ions  are  streaming  out  of  the  coupling  region 
faster  than  they  can  be  replaced.  This  distribution, 
containing,  a  'hole'  in  velocity,  is  known  to  suffer  an 
electrostatic  instability,  the  Post-Rosenbluth  'loss 
cone  instabi 1 ity. '  The  electrostatic  fields  which  are 
associated  with  this  instability  preferentially  pitch 
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angle  scatter  ions  with  high  transverse  velocities  (v  > 
vthermal)  into  the  parallel  direction,  where  they 
stream  out  the  loss  cone.  The  net  result  of  all  of  the 
above  is  the  establishment  of  a  quasi-steady  state  of 
equilibrium  consisting  of  a  continuous  flow,  out  of  the 
loss  cone  parallel  to  the  magnetic  field,  of  debris  and 
air  ions  with  velocities  significantly  larger  than  the 
transverse  thermal,  V ^  =  x  X  ~  2. -2. 5.  Thus  these 
escaping  ions  have  a  beam-like  character. “ 


At  a  working  group  meeting  at  MRC  in  September,  1979,  Palmadesso 
and  Dr.  Joseph  Workman  of  Berkeley  Research  Associates  suggested  a  method 
of  modelling  the  effects  of  the  Post-Rosenbluth  instability  in  our  SCATTER 
subroutine.  They  suggested  that  energy  be  removed  from  transverse  thermal 
at  the  rate 


d 

dt  ex 


(,  — \  vth 

-(ex  -  e)  _  if  ei  >  e, 

|  0 ,  i  f  ei  _<  ¥  . 


(3-6) 


where  ex  is  the  transverse  thermal  energy  density  at  a  particular  location 
in  the  blast  wave,  e  is  a  mean  thermal  energy  chosen  so  that  the  system 
will  approach  pressure  isotropization,  Vth  1s  the  local  mean  ion  thermal 
speed,  and  L  is  a  scattering  length,  proportional  to  the  Debye  length.  As 
this  energy  is  removed  from  transverse  thermal,  it  appears  as  parallel 
streaming  energy,  in  the  form  of  air  ions  traveling  at  parallel  velocities 
+  x  Vth.  with  X  =  2.5.  This  procedure  became  known  as  the  "Workman  - 
Palmadesso  formula,"  and  was  encoded  into  our  SCATTER  subroutine  for  run  B 
and  run  C.  The  effects  of  these  different  models  for  the  scattering 
process  are  explored  in  detail  in  Section  4. 
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The  other  change  which  has  been  made  in  CMHD,  and  which  was 
implemented  for  run  C,  was  prompted  by  informal  discussions  with  Dr. 
Palmadesso  of  NRL  and  Dr.  Steve  Brecht  of  SAI .  As  a  result  of  these 
discussions,  we  discovered  that  we  had  made  an  unfortunate  choice  of 
numerical  differencing  methods  for  an  equation  representing  the  "sling 
effect,"  mentioned  in  Section  2.  During  the  blast  wave  expansion  there  is 
interchange  between  energy  associated  with  the  expansion  transverse  to  $ 
(transverse  streaming  energy)  and  energy  of  ions  streaming  parallel  to 
The  relevant  differential  equations  are 


and 


^u,  ■  u,  .  V  6  .  S, 
^  v±  =  -uj  £  •  v  b*  .  I, 


(3-7) 

(3-8) 


where  U  is  the  parallel  ion  velocity,  V£  is  the  transverse  velocity  of 
of  the  plasma  (and  the  magnetic  field  lines)  and  \  and  e  are  unit  vectors 
parallel  and  transverse  to  B+,  respectively.  When  ions  streaming  parallel 
to  ti  find  themselves  on  moving,  curved  field  lines,  they  will  undergo  an 
interchange  of  energy  between  parallel  and  transverse  modes,  according  to 
these  equations.  This  is  simply  an  expression  of  the  fact  that,  under 
these  conditions,  the  unit  vectors  b  and  e*  are  changing  in  space  and  in 
time.  Total  energy  is  conserved,  of  course.  The  particular  numerical 
method  that  was  chosen  for  simulating  this  effect  had  the  unfortunate 
property  of  decreasing  the  parallel  energy  too  much,  under  certain  circum¬ 
stances.  After  this  discovery,  we  altered  the  code  to  correct  the  defi¬ 
ciency,  and  again  repeated  the  Starfish  calculation.  To  summarize,  then, 
there  have  been  three  CMHD  Starfish  calculations:  the  original  one, 
completed  in  1978,  and  referred  to  as  run  A,  which  utilized  the  diffusive 
form  of  the  SCATTER  subroutine;  run  B,  completed  in  1979,  which  used  the 
Workman  -  Palmadesso  (WP)  form  of  the  SCATTER  routine;  and  run  C,  which 
included  the  WP  SCATTER  routine  and  had  the  problem  with  the  "sling 
effect"  fixed.  The  three  calculations  used  identical  initial  and  boundary 
conditions. 
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Altitude  (km) 


Altitude  (km) 


The  overall  dynamics  of  the  blast  wave  expansion  were  affected 
very  little  by  these  changes  in  the  code.  Figure  5  and  6  show  comparisons 
of  electron  density  contours  at  0.3  seconds  and  at  1.2  seconds  after  burst 
time.  These  contour  plots  reveal  slight  differences  in  the  blast  wave 
expansions.  These  differences  are  almost  imperceptabie  in  the 
downward-going  section  of  the  blast  wave,  which  is  the  most  important  from 
the  standpoint  of  production  of  long-lasting  ionization.  The  differences 
in  blast  wave  expansion  are  more  obvious  in  the  horizontal  and  above¬ 
horizontal  directions,  but  are  well  within  the  range  of  uncertainty  which 
we  would  assign  to  any  one  of  the  calculations. 

The  motivation  for  performing  the  additional  CMHD  calculations 
was  to  determine  the  extent  to  which  the  debris  patch  predictions  would 
differ  from  that  of  the  original  Starfish  simulation.  Figure  7  shows  the 
time-integrated  energy  flux  predictions  for  the  downward  patch  as  a  func¬ 
tion  of  cross-field  radius.  Although  the  differences  are  well  within  the 
overall  range  of  uncertainty,  and  should  not  be  viewed  as  significant,  it 
is  worth  noting  that  the  three  calculations  differ  in  ways  that  are 
consistent  with  the  changes  that  were  made  in  the  code.  Runs  B  and  C  both 
result  in  higher  time-integrated  fluxes  at  small  radii  than  does  the 
original  run.  This  is  a  result  of  using  the  WP  SCATTER  subroutine,  which 
deflects  more  energy  into  high  parallel  velocities  during  the  early  stages 
of  the  blast  wave  expansion,  when  the  radius  of  the  blast  wave  is  small. 
The  curve  from  run  B  then  falls  below  that  for  run  A  at  larger  radii.  We 
might  expect  this  to  happen,  simply  because  if  more  energy  is  transported 
away  from  the  blast  wave  at  earlier  times,  there  is  less  available  at 
later  times  (and  larger  radii).  The  fact  that  curve  C  is  higher  than  the 
others  at  all  radii  is  a  consequence  of  the  changes  in  implementation  of 
equations  3-7  and  3-8,  which  resulted  in  keeping  more  energy  in  parallel 
modes  as  the  ions  moved  down  the  distorted  geomagnetic  field  lines. 
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Of  particular  interest  is  the  comparison  of  predictions  of  patch 
spectra.  Figures  8  through  10  show  velocity  spectra  of  energy  predicted 
to  arrive  in  the  downward  patch  within  100  km  of  the  burst  field  line. 
Energy  reaching  the  patch  in  the  form  of  fast  ions  is  plotted  in  Figure  8, 
fast  neutral  energy  is  shown  in  Figure  9,  and  the  sum  of  both  ion  and 
neutral  energy  is  given  in  Figure  10. 

The  two  more  recent  calculations  show  an  increase  in  the  amount 
of  energy  reaching  the  patch  at  high  velocities.  Run  C,  especially, 
predicts  substantially  more  energy  at  velocities  above  2000  km/sec.  In 
retrospect,  we  sucpect  that  the  changes  made  for  run  C  may  have  overcor¬ 
rected  the  earlier  deficiencies,  and  therefore  may  have  actually  caused  a 
bigger  boost  at  the  high  velocity  end  of  the  spectrum  than  is  warranted. 
The  qualitative  result  that  runs  B  and  C  predict  a  somewhat  harder  spec¬ 
trum  than  does  run  A,  is  to  be  expected,  in  view  of  the  changes  that  were 
made  in  the  CMHD  code.  The  low  velocity  end  of  the  spectrum,  however,  is 
essentially  the  same  for  all  three  runs.  This  is  an  important  result,  for 
two  reasons.  First,  it  demonstrates  that  the  exact  form  of  the  SCATTER 
subroutine  has  little  impact  on  this  portion  of  the  spectrum,  which  is 
comforting  in  view  of  the  fact  that  plasma  turbulence  is  essentially 
impossible  to  model  in  detail.  Also,  since  the  slower  ions  deposit  their 
energy  at  higher  altitudes,  where  the  resulting  ionization  will  tend  to 
persist  longer,  the  low  velocity  portion  of  the  patch  spectrum  is  more 
important  from  the  military  systems  viewpoint. 
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Figure  10.  Comparison  of  downward  (NCR)  debris  energy  patch  sped 
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SECTION  4 

PATCH  SPECTRUM  SENSITIVITY  STUDY 


To  obtain  more  detailed  information  about  the  sensitivity  of  the 
debris  energy  patch  spectrum  to  certain  features  of  the  physics  embodied 
in  CMHD,  we  performed  a  series  of  simplified  calculations  for  a  single 
magnetic  flux  tube.  We  chose  flux  tube  number  five  from  the  most  recent 
Starfish  simulation,  which  was  initially  vertical  and  included  the  volume 
between  cylindrical  radii  of  80  km  and  100  km.  The  flux  tube  geometry  and 
the  sources  of  energy  and  mass  were  taken  from  the  full  CMHD 
calculations. 

The  geometry  of  this  flux  tube  is  shown  in  Figure  11  at  four 
different  times  during  the  blast  wave  expansion.  At  .05  seconds  after 
burst  time,  the  flux  tube  has  been  compressed  and  distorted  in  the  hori¬ 
zontal  sector  and  in  two  or  three  sectors  above  and  below  horizontal.  At 
later  times,  more  of  the  flux  tube  is  involved  in  the  expansion,  and  by 
.80  seconds,  the  distortion  away  from  the  ambient  configuration  is  quite 
extreme,  involving  all  of  the  flux  tube  between  about  250  km  altitude  and 
1000  km  altitude.  As  the  blast  wave  compresses  the  flux  tube,  ions  and 
electrons  are  given  transverse  thermal  energy,  resulting  in  large  trans¬ 
verse  temperatures.  Collisional  processes  then  ionize  more  of  the  back¬ 
ground  air,  increasing  the  ion  mass  within  the  flux  tube.  The  rate  at 
which  transverse  thermal  energy  and  ion  mass  were  added  to  this  flux  tube 

were  determined  for  each  angular  sector  from  the  output  of  the  full  CMHD 

calculation.  These  rates,  integrated  over  all  sectors,  are  shown  in 

Figure  12.  The  time  integrals  of  these  rates  are  plotted  in  Figure  13. 
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igure  12.  Energy  and  mass  sources  for  the  entire  flux  tube,  as  dete 
mined  from  the  full  CMHO  run. 
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Figure  13.  Time  integrated  values  of  the  energy  and  mass  sources 


About  90%  of  the  energy  is  added  by  .10  seconds,  and  the  remaining  10%  is 
added  slowly  over  the  next  1.1  seconds.  There  is  a  peak  in  the  rate  of 
increase  of  ion  mass  at  .2  seconds,  but  it  is  much  less  pronounced  than 
the  peak  in  the  energy  source  function,  resulting  in  a  more  gradual 
increase  in  mass  over  the  running  time  of  the  problem.  It  is  worth  noting 
that  the  flux  tube  mass  increases  from  about  1x10s  gm  (ambient  ionosphere 
plus  ionization  produced  by  the  X-ray  flash)  to  about  60  times  this  value 
by  1.2  seconds.  As  we  will  see  later,  this  is  a  very  important  aspect  of 
the  phenomenology  of  this  burst. 

Having  obtained  the  flux  tube  geometry  and  the  sources  of  energy 
and  mass  from  the  output  of  the  full  CMHD  calculation,  we  then  did  a 
series  of  simplified  "CMHD-like"  calculations,  in  which  we  exercised  the 
SCATTER  subroutine  to  transfer  energy  from  transverse  thermal  into  paral¬ 
lel  streaming  energy,  and  the  TRANSPORT  subroutine  to  transport  ions  along 
the  distorted  field  lines  to  the  upper  and  lower  ends  of  the  flux  tube. 

The  block  diagram  of  Figure  14  illustrates  schematically  the 

flow  of  energy  and  mass  in  these  simplified  calculations.  As  energy  is 

added,  it  first  appears  as  transverse  thermal  energy.  Then,  provided  the 

turn-on  conditions  for  the  particular  "scattering  instability"  being 

modelled  are  met,  the  SCATTER  subroutine  illustrated  by  S  in  Figure  14 

shifts  ions  in  parallel  velocity  space,  from  small  and  large  U1  into 

larger  U  and  smaller  U  .  This  results  in  a  transfer  of  energy  from  E  to 

E  .  Ions  with  non-zero  U  then  stream  up  or  down  the  flux  tube.  This 
n  n 

latter  part  of  the  calculation  is  the  job  of  the  TRANSPORT  subroutine, 
illustrated  in  Figure  14  by  the  synbol  T.  In  a  particular  sector,  at  a 
particular  time,  the  amount  of  E|,  present  and  its  distribution  in  paral¬ 
lel  velocity  space  are  determined  by  the  time  integrated  effects  of  the 
scattering  process  in  that  sector  and  transport  between  that  sector  and 
those  above  and  below.  This,  in  turn,  determines  whether  or  not  the  turn¬ 
on  conditions  for  SCATTER  are  met. 
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igure  14.  Schematic  representation  of  the  flow  of  energy  and  mass  in  the 
single  flux  table  calculations.  Angular  sectors  5  through  21 
are  "active"  for  this  flux  tube  -  sectors  1  and  22  are 
straight  down  and  straight  up,  respectively,  where  the  non¬ 
orthogonality  of  the  CMHD  coordinate  system  causes  singulari¬ 
ties,  and  sectors  2  through  4  intersect  this  particular  flux 
tube  below  the  bottom  of  the  CMHD  grid,  ^represents  the 
action  of  the  SCATTER  subroutine,  which  takes  energy  from 
transverse  thermal  and  puts  it  into  paralle1  streaming,  and© 
represents  the  transport  of  parallel  strec  ng  energy  up  and 
down  the  flux  tube. 
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The  series  of  calculat ions  being  reported  here  was  done  to 
investigate  the  sensitivity  of  the  patch  spectrum  to  the  details  of  the 
particular  plasma  instability  model  encoded  in  the  SCATTER  subroutine. 

Results  from  eight  problems,  which  involved  variations  of  two  different 
instability  models,  are  shown  below. 

The  ions  that  reside  in  a  particular  spatial  location,  that  is, 
in  a  particular  angular  sector  of  a  given  flux  tube,  are  resolved  intc 
fifteen  groups  in  parallel  velocity  space,  covering  the  ranae  from  +4000 
km/sec  to  -4000  km/sec.  Velocities  parallel  to  are  considered  positive, 
antiparallel  velocities  are  negative  (the  undistorted  §  field  pointed 
downward,  as  at  the  geomagnetic  north  pole).  Initially  all  the  ions 
within  a  cell  are  in  the  U„  =  0  group,  as  shown  in  Figure  15.  As  the 

ion  mass  within  a  cell  increases  due  to  collisional  ionization  of  back¬ 

ground  air,  the  added  ions  also  appear  in  the  zero  velocity  group.  Then, 
if  the  instability  turn-on  condition  permits,  ions  are  "scattered"  from 
U„  =  0,  into  both  positive  and  negative  values  of  Uj ,  according  to  the 
particular  scatter  model  being  used.  The  two  types  of  models  used  here 
are  the  diffusive  scatter  model,  which  produces  a  distribution  function 

such  as  that  shown  in  Figure  16,  and  the  Workman-Palmadesso  (WP)  model, 
which  preferentially  scatters  high  velocity  ions  into  the  parallel  direc¬ 
tion,  as  is  illustrated  in  Figure  17. 

Figures  15  through  17  are  plots  of  ion  mass  per  unit  velocity  as 
a  function  of  parallel  velocity  for  one  angular  sector  within  a  flux 
tube.  By  combining  this  information  for  all  the  active  sectors,  we  can 
plot  contours  of  ion  mass  per  unit  velocity  as  a  function  of  parallel 
velocity  and  angular  sector,  resulting  in  displays  such  as  shown  in  Figure 
18.  The  detailed  results  of  the  present  calculations  can  best  be  summar¬ 
ized  by  a  sequence  of  plots  like  Figure  18,  and  a  series  of  these  plots  at 
succeeding  times  for  each  of  the  eight  calculations  will  be  shown  below. 


38 
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Figure  15.  An  example  of  the  initial  ion  distribution  function, 
within  one  angular  sector. 
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DMDU  VS  U  (SECTOR  11) 
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Figure  16.  An  example  of  the  distribution  function,  -ggp  within  one  angular 
sector,  after  the  diffusive  SCATTER  subroutine  has  deflected 
energy  from  the  transverse  direction  into  the  parallel  direction 
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DM  DU  (CM  /  (CM /SEC)  ) 


Figure  17.  An  example  of  the  distribution  function,  after  the  WP 

SCATTER  routine  has  deflected  energy  from  the  transverse  direction 
to  the  parallel  direction. 
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igure  18.  Typical  contours  of  the  logarithm  of  the  distribution  of  ion  mass 
per  unit  parallel  velocity,  with  dM/dUn  in  units  of  gm/(cm/sec), 
as  a  function  of  parallel  velocity  and  angular  sector  number.  The 
outermost  contour  has  dM/dUli  =  10”7  qm/(cm/sec),  and  the  innermost 
contour  has  dM/dUn  =  10  gm/(cm/sec~)  in  this  example  at  a  t ime  of 
0.1  sec. 


With  this  introduction,  then,  we  proceed  to  present  results  from 
the  eight  calculations.  The  first  calculation,  which  we  shall  call  CASE 
1,  used  the  diffusive  form  of  scattering  function,  which  performed  a 
simple  diffusion  in  velocity  space,  whenever  the  mean  parallel  thermal 
energy  per  ion  was  less  than  one-half  the  mean  transverse  thermal  energy 
per  ion.  That  is,  the  turn-on  condition  was 

e„  <  0.5  ei , 


with 


where 


i  °-5  WV2 

I  "V 

k 

1  mk  u.k 

l  mk 

k  K 


The  summations  are,  of  course,  over  the  15  parallel  velocity  groups.  The 
time  scale  for  the  diffusion,  or  scattering,  was  short  compared  to  other 
times  of  the  simulation,  so  that  the  value  of  e„  was  nearly  always  at 
least  as  large  as  it  would  be  for  complete  isotropizat ion. 


Figure  19  contains  plots  of  the  total  transverse  thermal  energy 
in  the  flux  tube,  the  total  parallel  energy,  and  the  energy  transported 
out  the  ends  of  the  flux  tube  as  functions  of  time.  The  transverse  energy 
peaks  at  about  .15  seconds,  thereafter  falling  gradually,  as  energy  is 
scattered  into  the  parallel  direction  at  a  faster  rate  than  it  is  being 
increased  by  continuing  compression  of  the  flux  tube.  The  parallel  energy 
rises  somewhat  less  rapidly  than  does  the  transverse,  peaking  at  about  .2 
seconds,  and  then  stays  fairly  constant.  This  points  out  an  important 
property  of  this  system,  namely,  that  since  the  plasma  instability  time 
scale  is  very  short  compared  to  other  time  scales  in  the  problem,  energy 
is  scattered  into  the  parallel  direction  at  whatever  rate  is  required  to 
maintain  a  quas i -equi 1 ibr i urn  between  the  scattering  process  and  transport 


of  parallel  energy  up  and  down  the  flux  tube.  That  is,  the  rate  of  trans¬ 
port  -  determined  by  the  geometry  of  the  flux  tube  and  the  distribution  of 
energy  in  parallel  velocity  space  -  is  the  controlling  rate.  The  third 
graph  in  Figure  19  gives  the  energy  which  has  been  transported  out  the  top 
and  bottom  of  the  grid,  and  subsequently  would  be  deposited  in  the  patch 
regions,  as  a  function  of  time.  It  rises  somewhat  more  steeply  at  first, 
when  the  fastest  ions  are  exiting  the  flux  tube,  and  somewhat  less  steeply 
later.  The  maximum  slope  occurs  at  about  0.3  seconds,  which  agrees  with 
the  time  of  the  peak  in  the  power-time  curve  inferred  from  optical  data. 

As  a  way  of  displaying  the  evolution  of  this  problem  in  more 
detail,  we  show  contour  plots  of  the  distribution  function  dM/dUn  at 
several  times  in  Figures  20  and  21.  The  problem  was  initialized  at  .025 
seconds,  and  the  correponding  ion  mass  distribution  function  is  strongly 
peaked  at  =  0.  The  peak  is  highest  in  sector  5,  which  is  at  the  lowest 
altitude  and  therefore  contains  the  most  X-ray-produced  ionization.  The 
contour  values  are  a  factor  of  ten  apart,  and  the  contours  vfich  appear  at 
this  time  range  from  the  lxlO'7  gm/(cm/sec)  to  lxlO~3  gm/(cm/sec).  At 
subsequent  times,  as  mass  and  energy  are  added,  and  as  the  diffusive  SCAT¬ 
TER  routine  is  turned  on,  the  distribution  function  is  spread  out  in  11, 
space,  and  the  peak  value  increases.  At  .050  seconds,  the  plot  is  fairly 
symmetric  about  U„  =  0,  due  to  the  symmetry  of  the  SCATTER  algorithm, 
and  the  limited  amount  of  elapsed  time  since  most  of  the  scattering  has 
occurred.  By  the  next  time,  however,  ion  transport  has  moved  those  parti¬ 
cles  with  anti-parallel  (upward)  velocities  to  higher  angular  sectors,  and 
those  with  parallel  (downward)  velocities  to  lower  sectors,  producing  a 
skewing  of  the  contours.  By  .20  seconds,  the  skewing  has  become  more 
pronounced.  At  progressively  later  times,  the  skewed  nature  persists,  but 
the  contours  retract  toward  the  center  as  the  faster  particles  are  trans¬ 
ported  out  the  ends  of  the  flux  tube.  As  this  is  taking  place,  the  SCAT¬ 
TER  routine  continues  to  diffuse  ions  in  parallel  velocity  space  as  needed 
to  keep  e,  up  to  the  "turn-on"  value.  However,  as  mass  is  added  -  the 
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total  ion  mass  increases  60-fold  during  the  course  of  the  problem  -  the 

average  transverse  energy  per  ion  becomes  smaller,  lowering  the  necessary 

spread  of  ions  in  Un  space.  At  1  second,  the  mass  in  the  problem  has 

increased  so  that  the  contours  range  from  IxlO-7  pm/(cm/sec)  to  lxio-2 

gm/(cm/sec).  The  calculation  was  continued  to  a  problem  time  of  2 

seconds,  at  which  time  the  simplifying  assumption  was  made  that  any 

further  evolution  of  the  distribution  function  would  be  due  solely  to 

transport.  We  could  then  construct  the  time-integrated  patch  spectra  for 

both  the  upward  -  and  downward  -  going  ions.  Figure  22  shows  the  spectrum 

for  the  downward  patch,  which  has  a  shape  character istic  of  CMHD  results, 

approximately  proportional  to  U_1. 

n 

CASE  2  utilized  the  WP  SCATTER  algorithm,  which  is  quite  differ¬ 
ent  from  the  diffusive  one.  Rather  than  simply  spreading  the  distribution 
function  out  in  Uj  space,  it  deflects  ions  directly  out  to  parallel 
velocities  given  by  ±  xVith  ,  where  the  NRL  recommended  value  of  x  is  2.5, 
and  V1  .  is  the  mean  transverse  thermal  velocity  of  the  ions  in  that  spa¬ 
tial  location.  The  total  transverse,  parallel,  and  patch  energies  as 
functions  of  time  are  plotted  in  Figure  23.  Figures  24  and  25  given 
contours  of  the  distribution  function  as  it  evolves  for  this  case.  The 
effects  of  the  WP  SCATTER  algorithm  are  quite  obvious,  particularly  at 
times  of  .05,  .10,  and  .20  seconds.  The  scattering  of  ions  out  to  L),  = 
x  V[th  ,  coupled  with  a  V[ that  increases  with  altitude  gives  a  "winged” 
appearance  to  the  contour  plots  at  .05  and  .10  seconds.  At  the  later 
times,  after  the  faster  ions  have  had  time  to  be  transported  out  of  the 
flux  tube,  these  effects  are  less  obvious.  The  resulting  patch  spectrum 
is  shown  in  Figure  26. 

CASE  3  used  a  combined  SCATTFR  algorithm,  that  both  diffused 
ions  in  velocity  space  and  deflected  ions  directly  out  to  +  x  V l  The 
other  five  cases  used  variations  of  the  WP  algorithm.  Table  I  lists  the 
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Figure  22.  The  velocity  spectrum,  -1^-,  for  the  downward  patch  for  CASE  1. 


Figure  24.  Distribution  function  contours  for  CASE 
and  .100  seconds. 
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Figure  26.  The  velocity  spectrum,  for  the  downward  patch  for  CASE  2. 
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Table  1.  Parameters  used  tor  the  eight  different  single  flux  tube  calcu- 
1  at  ions 


Type 

Value  of  x 

Turn-on 

Cond i t ions 

f 

mass 

Case  1 

Dif fus  ive 

M .  A . 

e  <  .5  e 

ii 

1.0 

2 

WP 

2.5 

e  -  .5  e 

l| 

1 .0 

3 

WP  +  jiff  us  ive 

2.5 

e,i  ■  .5  e 

1 .0 

4 

WP 

2.5 

e  <  10  e 

1.0 

5 

WP 

5.0 

e  .5  e 

ii 

e,i  '  -5  e 

1.0 

6 

WP 

2.5 

.10 

7 

WP 

2.5 

e  x  .5  e 

.01 

8 

WP 

2.5 

e„  "  -5  e 

0 

eight  runs,  the  type  of  scatter  algorithm  used,  the  value  of  x,  the  turn¬ 
on  condition,  and  the  value  of  fmass  ,  a  parameter  by  which  the  rate  of 
addition  of  ion  mass  was  multiplied. 

Comparisons  of  the  results  of  CASE  2,  CASE  3,  CASE  4  and  CASE  5 
are  made  in  Figures  27  through  33.  CASE  2  is  the  nominal  WP  scatter  run, 
while  CASE  3  includes  some  velocity  diffusion.  CASE  4  differed  from  CASE 
2  in  that  the  turn-on  condition  was  altered  to  allow  almost  all  (~  90%) 
of  the  available  thermal  energy  to  be  deflected  into  parallel  energy. 
CASE  5  was  different  from  the  nominal  case  in  that  a  larger  value  for  x 
was  used,  thus  deflecting  ions  to  higher  parallel  velocities  for  a  given 
V  Time  histories  of  transverse  thermal  energy  (E  )  for  the  four  cases 

are  plotted  in  Figure  27.  The  CASE  2  and  CASE  3  plots  are  very  similar. 
The  plot  for  CASE  4  shows  that  relatively  little  energy  was  left  in  E  , 
since  most  of  it  went  into  parallel  energy.  CASE  5  has  less  E  than  the 
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sverse  energy  as  a  function  of  time  for  CASE  2,  CASE  3,  CAS 


A 


first  two,  because  deflecting  ions  to  higher  parallel  velocities  results 
in  energy  being  transported  out  of  the  flux  tube  at  a  faster  rate,  thereby 
requiring  a  higher  rate  of  conversion  of  energy  from  to  Efl . 

Parallel  energies  are  plotted  for  the  four  cases  in  Figure  28. 
Again,  CASE  2  and  CASE  3  are  very  similar.  CASE  4  shows  substant  ial  ly 
more  E,  at  first,  but  after  about  .5  seconds,  it  has  about  the  same 
remaining  parallel  energy  as  the  others.  CASE  5  has  slightly  more  Efl  than 
the  first  two  cases  before  about  .3  seconds,  and  slightly  less  thereafter. 
This  difference  is  not  surprising,  since  CASE  5  deflected  ions  to  higher 
values  of  U ,  especially  early  in  the  problem,  which  subsequently 
increased  the  rate  of  transport  out  of  the  grid. 

Figure  29  shows  plots  of  patch  energy  vs.  time.  Again,  the 
first  two  cases  look  very  much  alike.  Cases  4  and  5  show  a  faster  rate  of 
delivery  of  energy  to  the  patch  at  early  times,  which  is  consistent  with 
the  above  discussion. 

Contours  of  the  distribution  functions  for  the  four  cases  at  .05 

seconds  are  presented  in  Figure  30.  Here,  we  see  a  difference  between 

cases  2  and  3.  The  latter,  having  employed  a  mixture  of  the  diffusive  and 

WP  scattering  algorithms,  shows  a  broader  distribution  near  U  =  0  than 

any  of  the  others.  CASE  4,  which  put  essentially  all  of  its  energy  into 

E^ ,  clearly  shows  more  ions  in  the  1000  km/sec  to  3000  km/sec  range  than 

the  other  cases  show.  CASE  5,  which  deflected  ions  out  to  U  =  ±  5  V  , 

i  i  tn 

shows  more  ions  at  the  extremes  of  the  velocity  range.  By  .4  seconds,  the 
time  of  the  distribution  function  plots  shown  in  Figure  31,  there  remain 
some  differences  between  the  various  cases,  but  they  are  much  less  obvious 
than  earlier.  Figure  32  shows  the  1  second  comparison.  At  this  time,  the 
differences  between  the  four  cases  are  insignif icant.  We  attribute  the 
similarity  in  these  plots  to  the  fact  that  the  energy  and  mass  source 
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Figure  28.  Parallel  energy  as  a  function  of  time  for  CASE  2,  CASE 
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functions  were  the  same  for  all  cases,  and  to  the  fact  that  the  control¬ 
ling  rate  is  that  of  transport,  rather  than  the  rate  of  scattering  in 
velocity  space.  Thus  the  details  of  the  scattering  mechanism  are  not  as 
critical  as  one  might  initially  presume  for  times  after  0.5  seconds  when 
the  low  velocity  patch  ions  are  generated. 

Figure  33  gives  the  patch  spectra.  The  most  significant  differ¬ 
ences  occur  at  velocities  above  3000  km/sec.  Clearly,  the  amount  of 
energy  which  reaches  the  conjugate  regions  at  these  high  velocities 
depends  on  the  plasma  instability  model  used  to  deflect  ions  into  the 
parallel  direction.  There  are  also  differences  in  the  low  velocity  (say, 
<  1500  km/sec)  portion  of  the  spectrum,  but  they  are  less  pronounced  than 
at  high  velocities.  The  one  thing  these  results  have  in  common  is  that 
they  all  show  a  substantial  low  velocity  component,  and  there  is  a  defin¬ 
ite  trend  for  the  energy  per  unit  velocity  to  increase  as  the  velocity 
decreases,  for  ll  <  2000  km/sec. 

tl 

The  final  comparison  we  shall  make  is  between  cases  2,  6,  7,  and 
8.  These  runs  differ  in  the  amount  of  ion  mass  added  during  the  running 
of  the  problem.  Of  course,  we  are  sure  that  CASE  2  has  the  appropriate 
amount  of  added  ion  mass,  since  we  are  confident  that  collisional  ioniza¬ 
tion  of  background  air  really  does  occur,  and  that  the  cross-sections  used 
are  reasonably  accurate.  The  purpose  of  this  comparison  is  to  demonstrate 
the  importance  of  this  aspect  of  the  early  time  phenomenology. 


Figures  34  through  36  show  comparisons  of  E(  ,  E  ( ,  and  Epatch- 
The  shapes  of  the  curves  for  the  four  cases  are  all  quite  similar.  The 


differences  indicate  that  as  the  rate  at  which  mass  is  added  is  reduced, 
the  rate  of  transport  of  energy  out  of  the  grid  is  increased,  thus  leaving 


progressively  less  energy  in  both  E;  and  E^for  cases  6,  7,  and  8  at  com¬ 
parable  times.  In  Figures  37  through  39,  we  compare  the  distribution 
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Comparison  of  distribution  function  contours  for  CASE  2,  CASE  6,  CASE 
and  CASE  8  at  1.00  seconds. 


functions.  At  .05  seconds,  the  differences  are  relatively  minor,  but  at 
.4  seconds  and  at  1  second,  there  are  substantial  differences.  The  most 
obvious  differences  are  that  the  reduced  mass  cases  show  contours  at 
larger  velocities  than  does  the  nominal  WP  CASE  2.  Not  quite  so  obvious, 
although  equally  important,  is  the  fact  that  as  the  added  mass  is  reduced, 
the  low  velocity  portion  of  the  distribution  is  reduced.  This  is  clearly 
evident  in  the  comparison  of  patch  spectra,  given  in  Figure  40.  The  mass 
added  to  the  flux  tube  through  collisional  ionization  of  background  air 
has  a  dramatic  effect  on  the  character  of  the  patch  spectrum.  Without 
this  added  ion  mass,  most  of  the  energy  arrives  in  the  patch  at  large 
velocities,  but  with  it,  the  patch  spectrum  peaks  at  low  velocities. 

Thus,  we  conclude  from  this  parameter  study  that  while  the  high 
velocity  portion  of  the  patch  spectrum  shows  substantial  dependence  on  the 
details  of  the  plasma  instability  model,  the  low  velocity  portion  is  rela¬ 
tively  insensitive  to  it.  On  the  other  hand,  the  entire  spectrum  is  sen¬ 
sitive  to  the  rate  at  which  background  air  is  ionized  in  the  expanding 
blast  wave  and  subsequently  shares  the  available  thermal  energy.  This  is 
a  fortunate  result,  because  we  have  a  great  deal  more  confidence  in  our 
ability  to  treat  collisional  ionization  properly  than  in  our  ability  to 
accurately  model  plasma  instabilities. 
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Defense  Nuclear  Agency 
ATTN:  STNA 

ATTN:  RAAE,  W.  McKechney 

ATTN:  RAAE,  H.  Fitz  Jr 

ATTN:  RAAE,  P.  Lunn 

3  cy  ATTN:  RAAE 

4  cy  ATTN:  TITL 

Defense  Technical  Information  Center 
12  cy  ATTN:  DD 

Dep  Under  Secretary  of  Defense 
Comm,  Cmd,  Cont  8  Intel  1 

ATTN:  Dir  of  Intelligence  Sys 

Field  Command 
DNA  Det  1 

Lawrence  Livermore  Lab 
ATTN:  FC-1 

Field  Command 
Defense  Nuclear  Agency 


ATTN: 

FCTT,  G.  Ganong 

ATTN: 

FCPR 

ATTN: 

FCTT,  W.  Summa 

ATTN: 

FCTT 

ATTN: 

FCTXE 

Interservice  Nuclear  Weapons  School 
ATTN:  TTV 

Joint  Chiefs  of  Staff 

ATTN:  C3S  Evaluation  Office  (HDOO) 

Joint  Strat  Tgt  Planning  Staff 
ATTN:  JLTW-2 


DEPARTMENT  OF  THE  ARMY 

Atmospheric  Sciences  Laboratory 
ATTN:  DELAS-EO,  F.  Niles 

BMD  Advanced  Technology  Ctr 

ATTN:  ATC-R,  W.  Dickinson 
ATTN:  ATC-R,  D.  Russ 
ATTN:  ATC-0,  W.  Davies 
ATTN:  ATC-T,  M.  Capps 

BMD  Systems  Command 

ATTN:  BMDSC-HLE ,  R.  Webb 
2  cy  ATTN:  BMDSC-HW 

Dep  Ch  of  Staff  for  Rsch  Dev  8  Acq 
ATTN:  DAMA-CSS-N 
ATTN:  DAMA-CSZ-C 
ATTN:  DAMA-WSZ-C 

Electronics  Tech  A  Devices  Lab 

ATTN:  DELET-R,  S.  Kronenberg 

Harry  Diamond  Laboratories 
2  cy  ATTN:  DELHD-NW-P 

US  Army  Ballistic  Research  Labs 
ATTN:  DRDAR-TSB-S 
ATTN:  DRDAR-BLB,  J.  Mester 

US  Army  Chemical  School 
ATTN:  ATZN-CM-CS 

US  Army  Corm-ElcC  Engrg  Instal  Agency 
ATTN:  CCC-EMEO-PED,  G.  Lane 

US  Army  Foreign  Science  S  Tech  Ctr 
ATTN:  DRXST-SD 

US  Army  Materiel  Dev  &  Readiness  Cmd 
ATTN:  DRCLDC ,  J.  Bender 

US  Army  Nuclear  8.  Chem  Agency 
ATTN:  Library 

US  Army  Research  Office 
ATTN:  R.  Mace 

US  Army  Satellite  Coitm  Agency 
ATTN:  Document  Control 

US  Army  TRADOC  Sys  Analysis  Actvy 
ATTN:  ATAA-PL 

USA  Missile  Command 

ATTN:  DRSM1-YS0,  J.  Gamble 


National  Security  Agency 

ATTN:  R-52,  J.  Skillman 

Under  Secretary  of  Def  for  Rsch  8  Engrg 

ATTN:  Strat  8  Thtr  Nuc  Forces,  B.  Stephan 
ATTN:  Strat  8  Space  Sys,  OS 
ATTN:  Defensive  Systems 

WWMCCS  Systems  Engineering  Org 
ATTN:  R.  Crawford 


DEPARTMENT  OF  THE  NAVY 

Naval  Electronic  Systems  Command 

ATTN:  PME  106-13,  T.  Griffin 
ATTN:  PME -117-2013,  G.  Burnhart 
ATTN:  PME  117-211,  B.  Kruger 
ATTN:  PME  106-4,  S.  Kearney 
ATTN:  PME  117-20 
ATTN:  Code  501A 


DEPARTMENT  OFJjH  E_N  AV  Y Jeon t. i nued ) 

Naval  Intelligence  Support  Ctr 
ATTN:  Document  Control 
ATTN:  N I  SC- 50 

Naval  Ocean  Systems  Center 

ATTN:  Code  5322,  M.  Paulson 
ATTN:  Code  5323,  J.  Ferguson 
ATTN:  Code  532 

Naval  Postgraduate  School 

ATTN:  Code  1424  Library 

Naval  Research  Laboratory 
ATTN:  Code  4187 

ATTN:  Code  6780,  J.  Fedder 

ATTN:  Code  4139,  D.  McNutt 

ATTN:  Code  4780,  S.  Ossakow 

ATTN:  Code  4780,  D.  Strobel 

ATTN:  Code  4720,  J.  Davis 

ATTN:  Code  6700,  T.  Coffey 

ATTN:  Code  7600,  B.  Wald 

ATTN:  Code  4700.1,  U.  All 
ATTN:  Code  2627 
ATTN:  Code  2000,  J.  Brown 
ATTN:  Code  6700 
ATTN:  Code  4128.2,  J.  Hohnson 
ATTN:  Code  7950,  J.  Goodman 

Naval  Surface  Weapons  Center 
ATTN:  Code  X211 

ATTN:  Code  F31 

Ofc  of  the  Dep  Ch  of  Naval  Ops 
ATTN:  OP  941D 

Strat  Systems  Project  Office 

ATTN:  NSP-2722,  F.  Wimberly 

ATTN;  NSP-2141 

ATTN:  NSP-43,  Tech  Lib 

DEPARTMENT.  OF  THE.  MR  FORCE 

Aeronautical  Systems  Division 

ATTN:  ASD/YH-YEF,  Capt  Guice 

Air  Force  Geophysics  Lab 
ATTN:  OPR- I 
ATTN:  OPR,  H.  Gardiner 
A1TN :  R.  O'Neil 
ATTN:  PHP 

ATTN:  PHI,  J.  Buchau 
ATTN:  R.  Babcock 
2  cy  ATTN:  OPR,  R.  Murphy 
2  cy  ATTN:  LKB,  K.  Champion 
2  cy  ATTN:  LKD,  R.  Narcisi 
2  cy  ATTN:  LKO,  R.  Huffman 
2  cy  ATTN:  OPR,  J.  Kennealy 
4  cy  ATTN:  CA,  A.  Stair 

Air  Force  Office  of  Scientific  Rsch 
ATTN:  AFOSR/NC 

Air  force  Systems  Command 


ATTN: 

DLAE 

ATTN: 

DLS 

ATTN: 

DLTW 

ATTN: 

DLXP 

ATTN: 

SDR 

DEPARTMENT  OF  THE  AIR  FORCE  (Co n t i nued j 
Air  Force  Tech  Applications  Ctr 


ATTN: 

TN 

ATTN: 

TT 

ATTN: 

TO 

ATTN : 

Tech  Library 

Air  Force  Weapons  Lab 
ATTN:  NTYC 
ATTN:  SUL 

Air  Force  Wright  Aeronautical  Lab 
ATTN:  A.  Johnson 
ATTN:  W.  Hunt 

Air  University  Library 
ATTN:  AUL-LSE 

Air  Weather  Service,  MAC 

ATTN:  DNXP,  R.  Prochaska 

Assistant  Chief  of  Staff 

Studies  S  Analyses 

ATTN:  AF/SASC,  C.  Rightmeyer 

Ballistic  Missile  Office 

AT1N:  ENSN  W.  Wilson 
ATTN:  SYC,  Col  Kwan 
ATTN:  SYDT 

Deputy  Chief  of  Staff 

Research,  Development  &  Aco 

3  cy  ATTN:  AFROS,  Space  Sys  &  C3  Dir 

Foreign  Technology  Division 
ATTN:  Nils  Library 

Rome  Air  Development  Cntr 
ATTN:  OCD,  J.  Simons 
ATTN:  TSLD 

Space  Division 
ATTN:  WE 

a’TN:  YGJB,  W.  Mercer 
ATTN:  YGP 

ATTN:  YKM,  Capt  Norton 
ATTN.  YN,  P.  Sivgals 

Strategic  Air  Command 

ATTN:  ADWAlf ,  R.  Bauer 
ATTN:  XPFS 
ATTN:  NRT 

DEPARTMENT  OF  ENERGY 

Department  of  Enerqy 
ATTN:  OMA 

OTHER  GOVERNMENT  AGENCIES 

Central  Intel  1 iqence  Agency 

ATTN:  OSWR/SSO  for  K.  Feuerpfetl 
ATTN:  OSWR/NLD 

Department  of  Commerce 

National  Bureau  of  Standards 
A1TN:  A.  Phelps 
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OTHER  GOVERNMENT  AGENCIES  (Continued) 


Department  of  Commerce 
National  Bureau  of  Standards 


ATTN: 

Sec  Ofc  for  M.  Krauss 

ATTN: 

Sec  Ofc  for  S.  Abramowitz 

ATTN: 

Sec  Ofc  for  J.  Devoe 

ATTN: 

Sec  Ofc  for  R.  Levine 

Department 

of  Commerce 

National  Oceanic  &  Atmospheric  Admin 

3  cy 

ATTN: 

E.  Ferguson 

3  cy 

ATTN: 

F.  Fehsenfeld 

Institute  for  Telecomm  Sciences 

ATTN: 

G.  Falcon 

ATTN: 

W.  Utlaut 

NASA 

ATTN 

J,  Watts 

ATTN 

W.  Roberts 

ATTN 

W.  Oran 

ATTN 

C.  Balcher 

ATTN 

N.  Stone 

NASA 

ATTN 

Technical  Library 

ATTN 

Code  6801,  A.  Tempkin 

ATTN 

Code  900,  J.  Siry 

3  cy 

ATTN 

A.  Aiken 

NASA 

ATTN 

P.  Kurzhals 

ATTN 

J.  Haughey 

ATTN 

A.  Schardt 

ATTN 

ST-5,  D.  Cauffman 

NASA 

ATTN 

MS- 168,  C.  Schexnayder 

NASA 

ATTN 

N-245-3,  R.  Whitten 

NASA 

ATTN 

J.  Gray 

DEPARTMENT 

0F_ENERGY  CONTRACTORS 

EG&G, 

Inc 

ATTN 

D.  Wright 

ATTN 

J.  Colvin 

University  of  California 

Lawrence  Livermore  National  Lab 

ATTN 

L-48,  E.  Woodward 

ATTN 

L- 325 ,  G.  Haugan 

ATTN 

L- 10,  H.  Kruger 

ATTN 

L-71,  J.  Chang 

ATTN 

L-31,  R.  Hager 

ATTN 

L-389,  R.  Ott 

ATTN 

Technical  Info  Dept  Library 

ATTN 

L-262,  W.  Duewer 

ATTN 

L -262 ,  D.  Wuebbles 

Sandia  National  Lab 

ATTN 

:  Org  12S0,  W.  Brown 

ATTN 

:  L.  Anderson 

ATTN 

:  Org  4231,  T.  Wright 

ATTN 

:  M.  kranm 

ATTN 

:  Tech  Lib  3141 

ATTN 

:  D.  Dahlqren 

DEPARTMENT  OF  ENERGY  CONTRACTORS  (Continued) 

Sandia  National  Labs,  Livermore 
ATTN:  T.  Cook 

Los  Alamos  National  Lab 


ATTN 

Librarian 

ATTN 

MS 362  Library 

ATTN 

T .  Kunkle ,  ESS-5 

ATTN 

M.  Sandford 

ATTN 

MS  670,  J.  Hopkins 

ATTN 

G.  Davis 

ATTN 

MS  664,  J.  Zinn 

ATTN 

D.  Simons 

ATTN 

R.  Jeffries 

ATTN 

J.  Wolcott 

ATTN 

G.  Barrasch 

ATTN 

R.  Carlos 

ATTN 

H.  Argo 

ATTN 

M.  Pongratz 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS 

Aero-Chem  Research  Labs,  Inc 

ATTN 

A .  Fon  t i j  n 

Aerodyne  Research,  Inc 

ATTN 

M.  Camac 

ATTN 

J.  Wormhoudt 

ATTN 

C.  Kolb 

Aerospace  Corp 

ATTN 

V.  Josephson 

ATTN 

R.  Slaughter 

ATTN 

J.  Straus 

ATTN 

II.  Mayer 

ATTN 

T.  Widhoph 

ATTN 

N.  Cohen 

ATTN 

J.  Reinheimer 

ATTN 

I.  Garfunkel 

ATTN 

1).  Olsen 

AVCO  Everett  Research  Lab,  Inc 

ATTN 

Aton 

ATTN 

f.  Von  Rosenherq  Jr 

ATTN 

Tech  Library 

Battele  Memorial  Institute 

ATTN 

H.  Lamuth 

ATTN 

STOIAC 

ATTN 

R,  Thatcher 

Berkeley  Research  Associates.  Inc 

ATTN 

J.  Workman 

ATTN 

S.  Brecht 

ATTN 

C.  Prettie 

Boeinq  Aerospace  Co 

ATTN:  MS/F7-63.  D.  Clauson 

Boston  Colleqe 

ATTN:  W.  Greider 

ATTN:  t.  Hrqhlom 

Boston  Colleup 

2  cy  ATTN:  Chairman  Dept  of  Chemistry 

2  cy  ATTN:  Chairman  Dppt  of  Physics 

University  of  California  at  Riverside 
ATTN:  A.  Lloyd 

ATTN:  0.  Pitts  Jr 


75 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


University  of  California  at  San  Diego 

Harris  Corp 

ATTN 

H.  Sooker 

ATTN 

E.  Knick 

Calspan  Corp 

Horizons  Technology,  Inc 

ATTN 

W .  Wurster 

ATTN 

R.  Kruger 

ATTN 

M.  Dunn 

ATTN 

J.  Grace 

HSS,  Inc 

ATTN 

C.  Treanor 

ATTN 

0.  Hansen 

Charles  Stark  Draper  Lab,  Inc 

Institute  for  Defense  Analyses 

ATTN 

J.  Gilmore 

ATTN 

H.  Gates 

ATTN 

A.  Tetewski 

ATTN 

H.  Wolf hard 

ATTN 

0.  Cox 

ATTN 

E.  Bauer 

University  of  Colorado 

International  Tel  4  Telegraph 

ATTN 

G.  Lawrence-LASP 

ATTN 

Technical  Library 

ATTN 

C.  Lineberger-JILA 

JAYCOR 

Columbia  University 

ATTN 

J.  Sperling 

ATTN 

Security  Officer  for  H.  Foley 

Kaman  Sciences  Corp 

Computer  Sciences  Corp 

ATTN 

J.  Jordano 

ATTN 

F.  Eisenbarth 

kaman  Sciences  Corp 

Comsat  Labs 

ATTN 

T.  Stephens 

ATTN 

D.  Fang 

kaman  Tempo 

Concord  Sciences 

ATTN 

W.  Knapp 

ATTN:  E.  Sutton 

ATTN 

B.  Gambill 

ATTN 

J.  Devore 

Cornell  University 

ATTN 

K.  Schwartz 

ATTN:  M.  kelly 

ATTN 

W.  McNamara 

ATTN:  D.  Farley  Jr 

ATTN 

J.  Devore 

ATTN 

V.  Stull 

University  of  Denver 

5  cy  ATTN 

DAS  I AC 

ATTN 

Sec  Officer  for  D.  Murcray 

Lockheed  Missiles  &  Space  Co, 

University 

of  Denver 

ATTN 

D.  Divis 

ATTN 

B.  Van  Zyl 

Lockheed  Missiles  4  Space  Co, 

E-Systems, 

Inc 

ATTN 

J.  Perez 

ATTN 

R.  Berezdivin 

ATTN 

J.  Cladis 

ATTN 

M.  Walt 

Electrospace  Systems,  Inc 

ATTN 

J.  Reagan 

ATTN 

P.  Phillips 

ATTN 

J.  Kumer 

ATTN 

B.  McCormac 

Environmental  Rsch  Inst  of  Michigan 

ATTN 

R.  Sears 

ATTN 

IRIA  Library 

Lockheed  Missiles  4  Space  Co, 

ESL,  Inc 

ATTN 

C.  Old 

ATTN 

J.  Marshall 

ATTN 

D.  Churchill 

ATTN 

R.  Heckman 

ATTN 

R.  Ibaraki 

University  of  Lowel 1 

ATTN 

E.  Tsui 

ATTN 

G.  Best 

ATTN 

J.  Lehman 

MIT  Lincoln  Lab 

General  Electric  Co 

ATTN 

D.  Towle 

ATTN 

R.  Edsall 

ATTN 

P.  Zav itsanos 

MA/COM  Linkabit  Inc 

ATTN 

1.  Jacobs 

General  Electric  Co 

ATTN 

H.  Van  Trees 

ATTN 

F.  Reibert 

ATTN 

A.  Viterbi 

General  Research  Corp 

Magnavox  Govt  4  Indus  Elec  Co 

ATTN 

J.  Ise  Jr 

ATTN 

G.  White 

ATTN:  8.  Bennett 

Martin  Marietta  Corp 

Geo-Centers 

ATTN 

R.  Heffner 

ATTN 

E.  Marram 
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McDonnell  Douglas  Corp 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


ATTN 

ATTN 

ATTN 


Technical  Library  Services 
R.  Halprin 
W.  Olson 


Meteor  Communications  Corp 
ATTN:  R.  Leader 

Mission  Research  Corp 


ATTN 

R.  Bigoni 

ATTN 

R.  Bogusch 

ATTN 

S.  Gutsche 

ATTN 

P.  Fischer 

ATTN 

F.  Guigliano 

ATTN 

G.  McCartor 

ATTN 

D.  Archer 

ATTN 

M.  Scheibe 

ATTN 

C.  Lauer 

ATTN 

R.  Hendrick 

4  cy  ATTN 

R.  Kilb 

4  cy  ATTN 

F.  Fajen 

5  cy  ATTN 

Tech  Library 

Mitre  Corp 

ATTN:  MS  0104,  M. 

ATTN 

B.  Adams 

Mitre  Corp 

ATTN 

W.  Hall 

ATTN 

J.  Wheeler 

ATTN 

W.  Foster 

Pacific-Sierra  Research  Corp 

ATTN:  H.  Brode,  Chairman  SAGE 

Photometries,  Inc 

ATTN:  I.  Kofsky 

Physical  Dynamics,  Inc 
ATTN:  A.  Thompson 

Physical  Oymanics,  Inc 
ATTN:  E.  Fremouw 

Physical  Research,  Inc 

ATTN:  R.  Deliberis 

Physical  Science  Lab 

ATTN:  W.  Berning 

Physical  Sciences,  Inc 


R&D  Associates 

ATTN 

W.  Karzas 

ATTN 

R.  Turco 

ATTN 

R.  Lelevier 

ATTN 

R.  Lindgren 

ATTN 

H.  Ory 

ATTN 

F.  Gilmore 

ATTN 

M.  Gantsweg 

ATTN 

B.  Gabbard 

Rand  Corp 

ATTN 

E.  Bedrozian 

ATTN 

C.  Crain 

Science  Applications,  Inc 

ATTN 

D.  Hamlin 

ATTN 

D.  Sachs 

ATTN 

L.  Linson 

Science  Appl ications ,  Inc 

ATTN 

J.  Cockayne 

Space  Data 

Corp 

ATTN 

S.  Fisher 

Spectral  Sciences,  Inc 

ATTN 

F.  Bein 

SRI  International 

ATTN 

G.  Smith 

ATTN 

R.  Tsunoda 

ATTN 

J.  Vickrey 

ATTN 

M.  Baron 

ATTN 

C.  Rino 

ATTN 

D.  McDaniels 

ATTN 

w.  Jaye 

ATTN 

V.  Gonzales 

ATTN 

R.  Hake  Jr 

ATTN 

R.  Leadabrand 

ATTN 

W.  Chesnut 

ATTN 

J.  Casper 

SRI  International 

ATTN 

C.  Hulburt 

Stewart  Radiance  Lab 

ATTN 

J.  Ulwich 

ATTN 

ATTN 

ATTN 


G.  Caledonia 
K.  Wray 
R.  Taylor 


Physics  International  Co 

ATTN:  Technical  Library 

Univ  of  the  Commonwealth,  Pittsburgh 


Sylvania  Systems  Group 

ATTN:  R.  Steinhoff 

Sylvania  Systems  Group 

ATTN:  J.  Concordia 
ATTN:  I.  Kohl  berg 

Technology  International  Corp 
ATTN:  W.  Boquist 

Teledyne  Brown  Engineering 


ATTN 

ATTN 

ATTN 


F.  Kaufman 
W.  Fite 
M.  Biondi 


ATTN: 

J. 

Ford 

ATTN: 

MS 

-12  Tech  Library 

ATTN: 

F. 

Leopard 

ATTN: 

N. 

Passino 

Princeton  University 

ATTN:  Librarian 

R&D  Associates 

ATTN:  B.  Yoon 
ATTN:  J.  Rosengren 


TRW  Electronics  &  Defense  Sector 
ATTN:  R.  Plebuch 

Wayne  State  University 
ATTN:  R.  Kummler 
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State 

University 

Visidyne,  Inc 

ATTN 

Sec  Con  ofc 

for 

K.  Baker 

ATTN:  W. 

Reidy 

ATTN 

Sec  Con  ofc 

for 

D.  Baker 

ATTN:  T. 

Degges 

ATTN 

Sec  Con  ofc 

for 

C.  Wyatt 

ATTN:  C. 

Humphrey 

ATTN 

Sec  Con  ofc 

for 

A.  Steed 

ATTN:  0. 

Shepa  rd 

ATTN 

Sec  Con  ofc 

for 

D.  Burt 

ATTN:  J. 

Carpenter 

ATTN:  H. 

Smi  th 

'Jlont 


’/Jayne  State  University 
ATTN:  W.  kauppila 


California  Institute  of  Technology 
ATTN:  J.  Ajello 


